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^ , Using improved theoretical calculations of the decay form factors in the Light Cone-QCD 

^ [ sum rule approach, we investigate the decay rates, dilepton invariant mass spectra and the 

'k^ ! forward-backward (FB) asymmetry in the decays B —> {K,K*)i~^i~ (£^ = e^,/x^,r^) in the 

;h ' standard model (SM) and a number of popular variants of the supersymmetric (SUSY) models. 

■ - - ■ Theoretical precision on the differential decay rates and FB-asymmetry is estimated in these 

theories taking into account various parametric uncertainties. We show that existing data on 
B -^ Xs7 and the experimental upper limit on the branching ratio B{B -^ K*^'^fi^) provide 
interesting bounds on the coefficients of the underlying effective theory. We argue that the 
FB-asymmetry in B ^ K*i~^i~ constitutes a precision test of the SM and its measurement in 
forthcoming experiments may reveal new physics. In particular, the presently allowed large- 
tan /? solutions in SUGRA models, as well as more general flavor- violating SUSY models, yield 
FB-asymmetries which are characteristically different from the corresponding ones in the SM. 
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1 Introduction 

Theflavor-changing-neutral-current (FCNC) transitions 5 -^ {Xs,Xd)'~f a.ndB -^ {Xs,Xd)i^i~ 
with Xs{Xd) being hadrons with overall strangeness S* = i^S* = 0), provide potentially strin- 
gent tests of the SM in flavor physics. FCNC transitions are forbidden in the SM Lagrangian 
and are induced by the GIM-amplitudes |l| at the loop level, which makes their effective 
strengths small. In addition, these transitions may also be parametrically suppressed in the 
SM due to their dependence on the weak mixing angles of the quark-fiavor rotation matrix — 
the Cabibbo-Kobayashi-Maskawa matrix VQj^jyj p|. These two circumstances make the FCNC 
decays relatively rare and hence vulnerable to the presence of new physics. In the context of 
the SM, the potential interest in rare 5-decays is that they would provide a quantitative de- 
termination of the quark-flavor rotation matrix, in particular the matrix elements Vtb, Vtd and 
^s 0) @) S il- A beginning in that direction has already been made by the measurement of the 
branching ratio B{B -^ X,7) J^, g, yielding \VtsV;^,\ = 0.035 ±0.004 §], in agreement with the 
expectations based on the CKM-unitarity 0]. Since complementary information will also be 
available from the B^-B'^- and B^-B^-mixing induced mass differences AMg and AM^, respec- 
tively, and from a number of rare kaon decays |Tl| , the parameters of the CKM matrix, which 



are already fairly constrained in the SM ||I2|, |1^, |TJ], will be multiply determined. This will 
result either in precise determination of the SM parameters in the flavor sector, comparable 
to the precision of the electroweak parameters of the SM [0, or, more optimistically, in the 
discovery of new physics. Thus, FCNC processes are potentially effective tools in searching 
for new physics, with the super symmetric theories receiving special attention in this context 

[Q, 0, |T3, ig, ig, gg 0, ^, 



Inclusive decay rates and distributions are relatively robust theoretically, making them 
well-suited to search for new physics which may result in distortions of the SM-distributions. 
Concerning rare -B-decays, we recall that the shape of the photon energy spectra in the radiative 
decays B -^ {Xs,Xd)'j depends on the underlying physics. However, deviations from the SM- 
based normalized photon-energy distributions are expected only for the low-to-intermediate 
photon energies, where the individual contributions from the various operators in the underly- 
ing effective theory are comparable. Measuring the low-E^ spectrum is, however, a formidable 
task in the present experimental set-up. More promising from the point of view of observing 
new-physics-induced distortions in the distributions are the decays B — > {Xs,Xd)i'^i~ , which 
provide the possibility of measuring Dalitz-distributions in a number of variables, which in 
turn could be used to determine the coefficients of the effective vertices in the underlying the- 



ory |[T^. This program is somewhat handicapped by the fact that heavy quark expansion in 



l/mfc breaks down near the end-points of the spectra P5| , ^ , near the cc-threshold and in 
the resonant region. Thus, a certain amount of modeling is unavoidable for the complete phe- 
nomenological profile of the decays B -^ Xgi^i^ ■ A number of studies has been undertaken 



to assess the non-perturbative effects [^, ^, ^ ^, ^, ^ 0, allowing to define limited 
kinematic regions where the short-distance physics in the SM and alternative theories can be 
quantitatively studied. 

While the inclusive rare decays discussed above are theoretically cleaner than exclusive 
decays, which require additionally the knowledge of form factors, they are also more difficult 
to measure. Present best limits from the CLEO collaboration on B —>■ XsfJ^^fi" and B —>■ 
Xse~^e~ [B3| decays are typically an order of magnitude larger than the corresponding SM- 



based estimates pH] . Moreover, inclusive rare decays are a challange for experiments operating 
at hadron machines. However, it is encouraging that the FCNC exclusive semileptonic decays, 
in particular the B -^ {K, K*)fi'^fi~ modes, are accessible to a wider variety of experiments. 
As we will argue quantitatively in this paper, some of the present experimental bounds on these 
(and related e^e~ modes) |^, 35] are already quite stringent. With the advent of the Fermilab 
booster, HERA-B, experiments at the LHC, and also the ongoing experiments at CLEO and 
the B-factories, the decays of interest B —>■ [K, K*)i^i~ will be precisely measured. It is 
therefore worthwhile to return to a comparative study of these decays in the SM and some 
candidate theories of physics beyond the SM to ascertain if these modes could be meaningfully 
used for searches of beyond-the-SM physics. 

In the context of the SM, exclusive FCNC semileptonic -B-decays have been studied in a 
number of papers fS^, ^, ^, |3^, 0, ^ ^, ^ Q with varying degrees of theoretical rigor 



and emphasis. The main purpose of this paper is twofold: First, we would like to report on 
an improved calculation of the decay form factors using the technique of the Light cone-QCD 
sum rules (LCSR) |^, ^]. Early studies of exclusive B decays in the LCSR approach were 
restricted to contributions of leading twist and did not take into account radiative corrections 
(see Refs. [J7, ^| for a review and references to original publications). In the present paper. 



we use the results of [|^ for vector form factors, which include NLO radiative corrections 
and higher twist corrections up to twist 4 [^ ^. For B ^ K form factors we improve 
on the results obtained in [^ by including the twist 4 mass correction terms calculated in 
1^. Second, we apply this technology to the SM and some popular variants of the SUSY 
models to determine the phenomenological profiles of the decays B -^ {K, K*)i~^i~ in these 
scenarios. For the latter, we choose minimal- and non-minimal SUGRA models, minimal 
flavor violating supersymmetric (MFV) model, and a general flavor-violating supersymmetric 
framework, studied in the mass insertion approximation (MIA). While all these models have 
been studied quite extensively in the literature for the inclusive decays B -^ Xsj and B -^ 
Xsi+t 



16, 17, IS, 19, P 



23|, IJ, we are not aware of corresponding studies for the 



exclusive decays. We strive to fill this gap in this paper. 

With our goals clearly stated, we turn to the main issues in the inclusive and exclusive 
rare 5-decays. Using the language of effective theories and restricting ourselves to the SM 
and SUSY, the short- distance contributions in the decays B — > Xgj and B —>■ Xsi^i~, and 

Cg and 



the exclusive decays of interest to us, are determined by three coefficients, called Cy*^^, 



Cio Q, 0-Q Of these, {Cj'^^l — the modulus of the effective coefficient of the electromagnetic 
penguin operator — is bounded by the present experimental measurements of the B -^ Xs'y 
branching ratio 0, ^. Using the 95% C.L. upper and lower bounds from the updated CLEO 
measurements 0]: 

2.0 X 10"^ < B{B -^ Xs-f) < 4.5 x 10"^ , (1.1) 

one gets in the next-to-leading precision the bounds. 



0.28 < ICj'^i 



mB)\ < 0.41 . 



:i.2) 



The magnitude of Cy^ (itib) in the SM [Q is well within the CLEO bounds but there is no 



experimental information on the phase of Cj^ (mB)- It is imperative to determine this sign 

"'^In general, more operators are present in supersymmetric theories and we discuss their possible effects later 
in this paper. 



experimentally, as it is model-dependent. In particular, in SUGRA-type theories, both positive 
and negative- valued solutions for Cj^ {ttib) are allowed in different SUSY-parameter regions. 
Despite the present lack of direct information on the sign of C7°^(mB), the bound in 
Eq. ( |1.2| ) is quite stringent and effectively limits possible new-physics effects due to the in- 
herent correlations among the branching ratio B{B — > Xg'y) and other observable quantities 
such as the B^-B^ mixing, ex and the mass of the CP-even Higgs boson, irth- In particular, in 
the context of the mSUGRA-models, present data on B{B —^ Xg'-^) 0, || and lower bounds on 
rrih 1^ do not allow the effective coefficient Cy^ {rnB) to have a positive sign [^. However, 



relaxing the GUT mass constraints on the parameters of the scalar superpotential, large-tan /5 
solutions exist, which are compatible with all present experimental constraints and predict a 



range of ?ti/i- values which are beyond the reach of LEP experiments p2| . Interestingly, these 
large-tan /^-solutions in non-minimal SUGRA models do admit positive values for C^^ {ttib) 
which are compatible in magnitude with the CLEO bounds. In a more general SUSY frame- 
work, the allowed parameter space for flavor- violating transitions is much larger. Thus, in the 
MIA- approach |^, not only the sign of Cy*^ but also that of Cio may have either value. As 
different dilepton invariant mass regions in B ^ Xsi^i~, the coefficients Cy^ (ttT'e), C*9^ ('^b) 
and Cio are weighted differently, a detailed knowledge of the invariant mass distribution and 



the FB-asymmetry 28[, together with the decay rate B -^ Xg'y, is completely sufficient to 



determine these effective coefficients fl^l-l With obvious changes, these remarks apply to 



the exclusive decays B -^ {K, K*)£~^i~ as well with the proviso that form factor dependence 
introduces an additional uncertainty, which we estimate in this paper. A relatively stable 
quantity is the value of the dilepton invariant mass for which the FB-asymmetry becomes zero 
in the SM. This has been discussed in the context of a number of phenomenological models 
for the form factors [Q. We argue here that using the results of the large-energy expansion 
technique (called LEET) ||5^, the uncertainty in the zero of the FB-asymmetry in the decays 
B -^ K*i~^i~ due to the form factors can be shown to be minimal. This yields a strikingly 
simple relation between the coefficients Cy^ and Cg^ which we present in this paper. 

This paper is organized as follows: In section 2, we introduce the effective Hamiltonian 
formalism for semileptonic rare 5-decays. Section 3 contains the definitions and derivations of 
the form factors in the decays B —>■ {K, K*)i^i~ using the Light cone-QCD sum rule approach. 
In section 4, we display the decay distributions for the invariant dilepton mass spectra for 
B — >■ {K, K*)i^i~ and the FB-asymmetry for B — > K*i^i~ . Section 5 contains our numerical 
results for the branching ratios and the FB-asymmetry in the SM, including comparison with 
the available data. Comparative studies in a number of variants of the supersymmetric models 
are presented in section 6. Section 7 contains a brief summary and some concluding remarks. 



^Note that Cj" , Cg and Cio sue Wilson coefficients (numbers), but Cg° is a function of the dilepton 
invariant mass and encodes also the information from the long-distance contribution. We assume that new 
physics leaves the long-distance part largely intact. 



2 Effective Hamiltonian 



At the quark level, the rare semileptonic decay b — » si^i can be described in terms of the 
effective Hamiltonian obtained by integrating out the top quark and W^ bosons: 

/-< 10 

In this paper, we use the Wilson-coefficients Ci calculated in the naive dimensional regulariza- 



tion (NDR) scheme |^ 



The above Hamiltonian leads to the following free quark decay amplitude: 



M{h^ str) 



^ VIV,, {C,^^ [s^.Lh] 



+ Cio [s7^L6] 



^75^ 



-2rhbCj''^ 



s 



i^n 



(2.2) 



Here, L/i?= (It 75)/2,s = 
respectively. We put msjrrH) 



g2, q 



p_(_ + ]9_ where p± are the four-momenta of the leptons, 
0, but keep the leptons massive. The hat denotes normalization 



in terms of the 5-meson mass, m^, e.g. s = s/m%, rhh = mb/rriB- Here and in the remainder 
of this work we shall denote by mi, = mb{fi) the MS mass evaluated at a scale fi and by rrih^poie 
the pole mass of the 6-quark. To NLO in perturbation theory, they are related by: 



rribdj,) = rribpoie 




-In — ^, 
4 fi^ 



(2.3) 



Note that A4{b ^ si^i~), although a free quark decay amplitude, contains certain long- 
distance effects from the matrix elements of four-quark operators, {i~^i~s\Oi\b), 1 < « < 6, 
which usually are absorbed into a redefinition of the short-distance Wilson-coefficients. To 
be specific, we define, for exclusive decay^, the effective coefficient of the operator Og = 

eV(167r2)(s7^L6)(l7'^^) as 



C9'="(s) = C9 + r(s), (2.4) 

where Y{s) stands for the above-mentioned matrix elements of the four-quark operators. A 



perturbative calculation yields |^, p5| 

^pert(s) = g{rhc,s) (3Ci - 



C2 + 3C3 + C4 + 3^5 + ^6 



-^(1, s) (4 C3 + 4 C4 + 3 C5 + Ce) - -^(0, s) (C3 + 3 C,) 



+ g(3C3 + C4 + 3C5 + C6) 



(2.5) 



We work in leading logarithmic (LLog) approximation with the values of Ci given in Table |1|. 
Formulae can be seen in |Q. For the decays B -^ Xs£~^£~ (likewise, for B -^ {K,K*)i~^i~), 



■^For inclusive decays one has in addition to take into account the 0{as) virtual and bremsstrahlung cor- 
rections to the matrix element {£~^ £' s\Og\b) as calculated in [ pO[ . For exclusive decays, one can define an 
effective coefficient by including only the virtual corrections. We do not include any perturbative corrections to 
the partonic matrix elements. However, corresponding corrections are included in the nonperturbative matrix 
element over mesons. 



Ci 


C2 


C; 


C4 


c. 


Ce 


Cf 


C9 


Cio 


cw 


-0.248 


+1.107 


+0.011 


-0.026 


+0.007 


-0.031 


-0.313 


+4.344 


-4.669 


+0.362 



Table 1: Values of the SM Wilson coefficients used in the numerical calculations, corresponding 
to the central values of the parameters given in Table \^. Here, Cj^ = C7 — C5/3 — Cq, and for 
Cg we use the NDR scheme and C(o) = 3Ci + C2 + SC^ + C4 + 3C5 + Cq. 



and with s far below the cc threshold, perturbation theory, augmented by power corrections, 
is expected to yield a reliable estimate. The power corrections in 1/ml can not be calculated 
near the threshold s = 4m^ and in the resonance regions, as the heavy quark expansion 



breaks down [^. So, a complete profile of the FCNC semileptonic decays can not at present 
be calculated from first principles. Several phenomenological prescriptions for incorporating 
the nonperturbative contributions to Y{s) exist in the literature |2^ 



29| , |31| . The resulting 



uncertainties on Cg^ and various distributions in the inclusive decays have been worked out 



m 



i 



30 



to which we refer for detailed discussions. 

m 



In the present paper we use the 
and Ali, Mannel and Morozumi pB| , 



two parametrizations due to Kriiger and Sehgal 

and interpret the difference in results for Cq^^ as an estimate of the theoretical uncertainty. 
Nonperturbative effects originate in particular from resonance corrections to the perturbative 
quark-loops included in Fpert(s). Light-quark loops are suppressed by small Wilson-coefficients, 
so it is essentially only the charm-loop that matters. Ref. [^] suggests to add the cc resonance- 
contributions from J/\&, \E'', . . . , ^E'^^^ to the perturbative result, with the former parametrized 
in the form of a phenomenological Breit- Wigner Ansatz |^ . Y is then given by 



K, 



-' pert l^"'. 



+ ^cw 



E 

V,='4>{ls),... 



ri<i 



r(v. 



my. 



,i/.(6s) 



my/ 



smB^ — imy.T 



Vi 



(2.6) 



with C(°) = 3Ci + C2 + 3C3 + C4 + 3C5 + Ce. 
factorization approximation which with Nc 
too small branching fraction for B -^ K'^*'>Vi. 



The phenomenological factors /?» correct for the 
= 3 (also called naive factorization [^) gives a 
They can be fixed from 



B{B -^ K^*'>Vi)B{V, 



(2.7) 



where the right-hand side is given by data |T^. While in the literature for inclusive B -^ 
Xsi^i~ decays, one comes across a universal K.i{Xs) = ni{Xs) = 2.3 [^, we have evaluated 
the individual factors for the lowest two cc resonances, shown in Table ^ In our numerical 
analysis we use for the higher resonances \Ef(**), . . . , \If(^) the average of J/\E' and \E''. We have 
averaged over charged and neutral B mesons if data are available. Concentrating on J/"^ only 



K 


J/m 


m' 


K 
K* 


2.70 
1.65 


3.51 
2.36 



Table 2: Fudge factors in B 
form factors. 



K^*)j/m^m' -^ K'-*k+£- decays calculated using the LCSR 



and assuming that the inclusive case Xs is saturated by K and K* , we get ki(Xs) = 1.9. Note 
that only the combination |C*^''^/i:i| can be fixed from the J/?/', T/^'-data. However, we treat the 
phase of the Kj as fixed to the one in the factorization approach. This is substantiated by data 
in which the Bauer- Stech-Wirbel parameters ai and 02 [Q are consistently determined, with 
Oi coming out close to its perturbation theory value and the sign of 02/01 is the one given by 
the factorization approach [^, pT| . 

In the AMM-approach, it is tacitly assumed that the extrapolation of the Breit-Wigner 
form away from the resonances could be used to estimate these power corrections reliably. 

The KS-approach |2^, on the other hand, bypasses the perturbative/non-perturbative di- 



chotomy by using the measured cross-section a{e^e~ — >■ hadrons) together with the assumption 
of quark-hadron duality for large s to reconstruct Y{s) from its imaginary part by a dispersion 
relation. However, perturbative contributions in cr(e^e~ -^ hadrons) and B — >■ Xs£~^i~ are not 
identical. In particular, the perturbative part of Y{s) has genuine hard contributions propor- 
tional to ml, which can neither be ignored nor taken care of by the quark-duality argument. 
The issue in this approach remains as to how much of the genuine perturbative contribution in 
B -^ Xsi~^i~ arising from the cc-continuum should be kept and there is at present no unique 
solution to this problem, as argued in M^ to which we refer for further discussion of this 



point. As stated earlier, we shall take the difference between the AAM-based and KS-based 
approaches for the long-distance contributions as a theoretical systematic error. 

3 Form factors from QCD sum rules on the light-cone 

Exclusive decays B — »• [K, K*)i~^i~ are described in terms of matrix elements of the quark 
operators in Eq. ( p.2|) over meson states, which can be parametrized in terms of form factors. 
Let us first define the form factors of the transition involving the pseudoscalar mesons 
B ^ K. The non-vanishing matrix elements are {q = Pb — p) 

{2 2 ^ 22 
nfYl HI I HI HI 

iPB+p),- \ ^ g,|+ \ '' fo{s)q„ (3.1) 

and 

{K{p)\sa,,q%l + ^,)b\B{pB)) ^ {K{p)\sa,,q%\B{pn)) 

= ^{{PB+p)^,s-q^,iml-mj^)} . (3.2) 

•■ J rriB + rriK 

For the vector meson K* with polarization vector e^, we can define the semileptonic form 
factors of the V — A current by 

{K*{p)\{V-A)^\B{pb)) = -ie*JmB + mK*)A,{s)+t{pB+p),{e*PB) "^'^'^ 



rriB + rriK* 



+iq,{^*PB) ^^^^ {A,{s) - A,{s)) + e^.,.e*»'^ '^^}'^ . (3.3) 

s rriB + rriK* 



Note the exact relations 



mB + rriK* . , . ttib - tuk* . , . 

A(0) = ^3(0), 
{K*\d,A^^\B) = 2mK*{e*pB)Ao{s). (3.4) 



The second relation in (p.4| ) ensures that there is no kinematical singularity in the matrix 
element at s = 0. The decay B — > K*i~^i~ is described by the above semileptonic form factors 
and the following penguin form factors: 

{K*\sa^,q''{l + ^,)b\B{pB)) = ze^,,.e»^ 2Ti(s) 

+ T,{s) {e;(m| - ml,) - {e*pB) {pb + p),} 

+ T3(s)(e>B) W/. 2 —iPB+p)f,\ (3.5) 



m^ - rriK* 



with 



Ti(0) = r2(0). (3.6) 



All signs are defined in such a way as to render the form factors positive. The physical range 
in s extends from Smjn = to Smax = ("^b — ttT'K,k*)^- 

Lacking a complete solution of non-perturbative QCD, one has to rely on certain approxi- 
mate methods to calculate the above form factors. In this paper, we choose to calculate them 
by the QCD sum rules on the light-cone (LCSRs). The method of LCSRs was first suggested 
for the study of weak baryon decays in [^ and later extended to heavy meson decays in P6 |. 



It is a nonperturbative approach which combines ideas of QCD sum rules [^ with the twist 
expansion characteristic for hard exclusive processes in QCD [^ and makes explicit use of 
the large energy of the final state meson at small values of the momentum transfer to leptons 
s. In this respect, the LCSR approach is complementary to lattice calculations |E^, which 
are mainly restricted to form factors at small recoil (large values of s) and at present require 
the scaling behavior found in the context of the LCSRs to extrapolate to smaller values of s 
[ ^ . Of course, the LCSRs lack the rigor of the lattice approach. Nevertheless, they prove 
to provide a powerful nonperturbative model which is explicitly consistent with perturbative 
QCD and the heavy quark limit. 

Early studies of exclusive B decays in the LCSR approach were restricted to contributions 
of leading twist and did not take into account radiative corrections. These corrections, included 
in the estimates presented here, turn out shift the form factors by ~ 10%. 

In previous works |]69| , ^, ^, the resulting form factors have been parametrized by a 
modified single pole formula, 

1 — aps + bps^ 

obtained from a fit to the LCSR result in the region s < 0.54. The extrapolation of this 
parametrization to maximum s is prone to spurious singularities below the physical cut starting 





/+ 


/o /t 


Ai A2 Ao V 


Ti T2 T3 


F(0) 


0.319 


0.319 0.355 


0.337 0.282 0.471 0.457 


0.379 0.379 0.260 


Cl 


1.465 


0.633 1.478 


0.602 1.172 1.505 1.482 


1.519 0.517 1.129 


C2 


0.372 


-0.095 0.373 


0.258 0.567 0.710 1.015 


1.030 0.426 1.128 


C3 


0.782 


0.591 0.700 









Table 3: Central values of parameters for the parametrization (^.Tj) of the B ^ K and B -^ K* 
form factors. Renormalization scale for the penguin form factors fx and Ti is fi = mi,. C3 can 
be neglected for B — >■ K* form factors. 





/+ 


/o fr 


Ai A2 Ao V 


Ti T2 T3 


F(0) 


0.371 


0.371 0.423 


0.385 0.320 0.698 0.548 


0.437 0.437 0.295 


Cl 


1.412 


0.579 1.413 


0.557 1.083 1.945 1.462 


1.498 0.495 1.044 


C2 


0.261 


-0.240 0.247 


0.068 0.393 0.314 0.953 


0.976 0.402 1.378 


C3 


0.822 


0.774 0.742 









Table 4: Parameters for the maximum allowed form factors. 





/+ /o It 


Ai A2 Ao V 


Ti T2 T3 


F(0) 


0.278 0.278 0.300 


0.294 0.246 0.412 0.399 


0.334 0.334 0.234 


Cl 


1.568 0.740 1.600 


0.656 1.237 1.543 1.537 


1.575 0.562 1.230 


C2 


0.470 0.080 0.501 


0.456 0.822 0.954 1.123 


1.140 0.481 1.089 


C3 


0.885 0.425 0.796 









Table 5: Parameters for the minimum allowed form factors. 
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Figure 1: LCSR form factors with theoretical uncertainties for the B —^ K transition as a 
function of s. Solid, dotted and dashed curves correspond to /+, /r, /o, respectively. Renormal- 
ization scale for /r is fi = rub- 



(b) 





0.5 



Figure 2: LCSR form factors with theoretical uncertainties for the B -^ K* transition as a 
function of s. In (a), the solid, dotted, dashed and short long dashed curves correspond to 
V, Aq, Ai, A2 and in (h), the solid, dotted and dashed curves correspond to Ti,T2,T3, respec- 
tively. Renormalization scale for Ti is n = mi,. 



at s = m^. . In the present work we thus choose a different parametrization which avoids this 
problem: 

F{s) = F(0) exp(cis + css^ + cgs^^). (3.7) 

The term in s^ turns out to be important in i? ^ i^ transitions, where s can be as large as 
0.82, but can be neglected for B -^ K* with s < 0.69. The parametrization formula works 
within 1% accuracy for s < 15GeV^. For an estimate of the theoretical uncertainty of these 
form factors, we have varied the input parameters of the LCSRs, i.e. the b quark mass, the 
Gegenbauer-moments of the K and K* distribution amplitudes and the LCSR-specific Borel- 
parameters M^ and continuum threshold sq within their respective allowed ranges specified in 
[0, ^ and obtain the three sets of form factors given in Tabs. 0-^, which represent, for each 



s, the central value, maximum and minimum allowed form factor, respectively. We plot the 
form factors in Figs. |l| and ^ 

Our value of Ti(0) is consistent with the CLEO measurement of B{B -^ K*'y)exp = (4.2 ± 
0.8 ± 0.6) • 10^^ ||70[. From the formula for the decay rate. 



T{B -. K*^) = ^i^§^mlmUl - mlJmir\Cr^^mm' , (3-8) 

the central values of the parameters given in Table §, Ti(0) = 0.379 and with tb = 1.61 ps we 
find B{B -^ K*-f)th = ^A ■ 10"^. 

4 Decay Distributions 

In this section we define various decay distributions whose phenomenological analysis will be 
performed in the next section. 
Eq. ( |2.2| ) can be written as 

M = ^y:syt,ms [r; {li^i)+r'^ {iri.i)] , (4.1) 

where for B — > Ki~^i~, 

T; = A'is)p, + B\s)q,, (4.2) 

T^' = C'{s)p, + D'{s)q,, (4.3) 

and for B -^ K*i+i-, 

T; = A{s) e^,,^e*^p^4, - iB{s) t\ + iC{s) (e* • pb)p, + iD{s) (e* • pb)% , (4.4) 

T^ = E{s) e,,^pe*'p%p^^. - zF(S) e% + tG{s) (e* ■ Pb)p, + tH{s) (e* • pB)q, , (4.5) 

with p = pb + Pk,k*- Note that, using the equation of motion for lepton fields, the terms in 
q^ in T^ vanish and those in T^ become suppressed by one power of the lepton mass. This 
effectively eliminates the photon pole in B' for B -^ K. 
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The auxiliary functions above are defined as 
A'is) = C^^is) Uis) + 



2^^ Cr^^Ms) 



1 + rhK 



1 + rriK* s 



B(s) = (l+mK*] 



C{s) 



1 - mj^. 



s 



(1 - mK*)C,''''is)A2{s) + 2m,C,'^ (t,{s) + L^^T2(s) 



D{s) = i [C^'^{s) ((1 + mK^)Ai{s) - (1 - mK')A2{s) - 2mK*Ao{s)) 



-2mbCr'''T,(s] 



-C,oV{s) 



E{S) = z . 

l+rriK* 

F{s) = {I + mK^)CioA,{s) , 

1 



G{s) 



1 +rhK* 



-CioA2(s) , 



(4.6) 

(4.7) 

(4.8) 
(4.9) 

(4.10) 

(4.11) 

,(4.12) 



(4.13) 
(4.14) 
(4.15) 
(4.16) 



H{s) = ^Cio [(1 + mxOA(s) - (1 - mx0^2(s) - 2m^. Ao(s)] . 



(4.17) 



Note that the inclusion of the full s-quark mass dependence into the above formulae can be 
done by substituting mi, — > m^ + rris into all terms proportional to C-j'^ Ti and Cj'^ fx and 
mh-^mh- iTLs in Cy'^^Ts^s, since 0^ ~ sa^^, [(m^ + mj + (mj, - 771^)75] q^h. 
We choose the kinematic variables is, u) to be 



which are bounded as 



s = q^ = {p++p^f , 

u = {pB-p-f - {pB-p+f 

(2fnif < s <{l-rnK,K')\ 
—u{s) < u < u{s) , 



(4.18) 
(4.19) 

(4.20) 
(4.21) 
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with rn£ = mi/niB and 



u\s\ 



'A(l-4 



nip 



X = A(l, rhK,K*^ s) = 1 + rfiK^K* + S^ - 2s - 2m^_^^(l + s) 



(4.22) 



(4.23) 



Note that the variable u corresponds to 6, the angle between the momentum of the i?-meson 
and the positively charged lepton i~^ in the dilepton CMS frame, through the relation u = 
—u{s) cos 6 p8[. Keeping the lepton mass, we find the double differential decay widths T^ and 
r-'^* for the decays B -^ Ki^i~ and B — > K*i~^i~, respectively, as 



d^r^ GWm 



dsdu 






2i%5 
x{(|Af + |Cf)(A-M2) 

+ |Cf4m/(2 + 2m|^ - s) + Re{C'D'*)8W{l - ml) + |Df 4m/s} , (4.24) 



d^r^* G|«2^5 



dsdu 



X- 



+ 



21%^ 



^wzv,,:^ 



'ts "^M 



E\ 



4rh|^. 



(s(A + u^) + Am^^X) + ^ (s(A + u^) - Am^^X) 



|5p (a - M^ + 8m^, (s + 2mj)) + jFj^ (a - m^ + 8m^, (s - Am 



- 2sM [Re(SE*) + Re{AF*)] 
X 



+ 



Am\, 

1 
2m\, 



\C\\X - u^) + |G|2 (a - m' + 4m2(2 + 2rn\, 
Re(EC*)(l - m^. - s)(A - u^) 



+Re{FG*) (1 - m|^* - s) (A - n') + 4m|A 



m; 



mt^. 



_ft:* 



Re(Fi7*)-Re(GiJ*)(l-m^*) + |/7|'-r^sA 



m; 



m 



(4.25) 



ic* 
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4.1 Dilepton mass spectrum 



We now give formulas for the dilepton invariant mass spectra. Integrating over u in the 
kinematical region given in Eq. ( |4.21|) we find 



dr 



K 



Gla^m% ,^^^^^ ,2, 



ds 



210^^5 



V:sVt,\u{s) 



x{(\A'\' + \C'\'){X 



U[S 



+ |Cf 4m/(2 + 2ml - s) + Re{C'D'*)8me\l - m^) + |Df 4m/s} , 



(4.26) 



dr 



K* 1^2 ^2 5 



X' 



"2 " / ^n2 



sA(l + 2-:^) + |E|"s- 



^fn\. 



\B\\X - ^ + 8m^.(5 + 2m2)) + |F|2(A - ^ + 8m|,,(S - 4m,2)) 



imj^. 



\C\\X- 



uis 



IGP A 



■us 



4m^(2 + 2m^. - s) 



2m|^* 



Re(5C*)(A 



uis] 



-){l-ml,~s) 



+ Re(FG*)((A 



Uis 



-){l-ml, -s) +4m^A) 



-2^A 



m 



K* 



Re{FH*) - Re{GH*){l - rh\,\ 



+ -^s\H 



m 



K* 



(4.27) 



Both distributions agree with the ones obtained in |^. In the limit m^ — i> the form factors 
/o (or /_) and Aq do not contribute. Furthermore, since \Cj'^^\ <ti \C<^'^^\, [C'lol, the influence 
of fxyTs on the distributions is subdominant. That means that roughly dV^/ds ~ |/+P for 
i = e, fi in the low s region below the J/\l/, with a ~ —12% effect coming from Cy^ fr terms. 
For B -^ K*, the b ^ S'j transition is more important: for s < IGeV^ the photon pole is the 
dominant contribution, and it still contributes ~ —30% around s ~ 3GeV^. 



4.2 Forward- backward- asymmetry 



The differential forward-backward-asymmetry (FBA) is defined as |l28 

d^FB 



ds 



d-uds 



d^r 

-u{s) duds 



(4.28) 
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The FBA vanishes in B ^ K£^i decays as can be seen from Eq. ( |4.24] ), since there is no 
term containing it with an odd power. For B — > K*i'^i~ decays it reads as follows 

d^FB Gpa m^ |t/*t/ |2 --/-N2rT3 /m?*\ I T? ( ^T?*\^ 
— \^ts^tb\ su{s) [tie{BL ) + tie{Al' )\ 



ds 2^^TT^ 

xCio 



^' (4.29) 



s 



The position of the zero Sq is given by 



cSr-^^ ^b^es\T2{so) ,, _._ , , Ti(so) 



Re(C9^"(so)) = -^C/" -^(l-m;,0 + -r^(l + ^i.O , (4-30) 

■So t^u-soj ^ [So) J 

which depends on the value of m^, the ratio of the effective coefficients Cj'^^ /Re{Cg'^^{so)), and 
the ratio of the form factors shown above. It is interesting to observe that in the Large Energy 
Effective Theory (LEET) ^8|, both ratios of the form factors appearing in Eq. ([4.30|) have 
essentially no hadronic uncertainty, i.e. , all dependence on the intrinsically nonperturbative 
quantities cancels, and one has simply: 

T2 1 + TTlK* ( S \ 



Ai 1 + rhr^t — s \ \ — rh 



Ti ^ 1 

V 1 + rhx* 



(4.31) 



With these relations, one has a particularly simple form for the equation determining sq, 
namely 

Re{C,^^{so)) = -2'^C,^ \J-'' . . (4.32) 

So 1 + mj^, - So 



Thus, the precision on the zero-point of the FB-asymmetry in i? ^ K*vi~ is determined 
essentially by the precision of the ratio of the effective coefficients and rrii,, making it at par 
with the corresponding quantity in the inclusive decays B — >■ Xsi^i~ , for which the zero-point 
is given by the solution of the equation Re(C9'^^(so)) = —^Cy'^^. We find the insensitivity of 
So to the decay form factors in i? ^ K*i~^i~ a remarkable result, which has also been discussed 
in [^. However, the LEET-based result in Eq. ( |4.31| ) stands theoretically on more rigorous 
grounds than the arguments based on scanning a number of form factor models. With the 
coefficients given in Table |1| and nib = 4.4 GeV, we find so = 0.10 (i.e. so = 2.9 GeV^) in the 
SM. From Eq. ( [4.30| ) it follows that there is no zero below the cc reonances if both Cg and Cj^ 
have the same sign as predicted in some beyond-the-SM models. 

From the experimental point of view the normalized FB-asymmetry is more useful, defined 

as 

d^FB _ d^FB ;dr 

ds ds ds 
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which is equivalent to the energy asymmetry |]TB|, ^. A shghtly different definition is 

d^PB d^FB 



ds 



ds 



-/r 



(4.34) 



whose integral gives the global energy asymmetry A'p^ = Afb/^- 
We summarize the characteristics of our observables: 



^(B -^ Ki+e , ^..^ , 

as ^ 'as 

sign(C7^^i?e(C9'=^)) = +1. 



and ^{B 



K*£'^i ) get maximal for maximal \C' 



cffi 



IC9I, |Cio| and 



J^'^ {B — > K*£~^£~) is proportional to Cio and has a characteristic zero (barring the 
trivial solution Cio = 0, which we do not entertain here) if Eq. (|4.3CI|) is satisfied, which 
requires 



cS\ 



signiCr'^^ReiC''^)) 



(4.35) 



The condition in Eq. ( ^.35| ) provides a discrimination between the SM and models having new 
physics. For example, this condition is satisfied in the SM and the SUGRA models with low- 
tan /3, in which case the actual position of Sq would provide the further discriminant. However, 
it turns out that the allowed parameter space of the SUGRA models with large-tan /5 yield 
sign(C7'^ Re{Cg'^ )) = +1 [^, leading to the result that the FB-asymmetry in these models is 
parametrically different. In particular, in all such cases, there is no zero of the FB-asymmetry. 

5 Branching Ratios and FB- Asymmetry in SM 

The input parameters that we use in our numerical analysis are given in Table ^ The param- 
eters which are either well-known or have a small influence on the decay rates have been fixed 
to their central values, but we vary four of the listed parameters, nit, ^, rub^poie and ^^(m^), 
in the indicated range. Furthermore, in the evaluation of the various distributions we use for 
rhb the MSbar mass evaluated at the scale /i = rrib^poie, see Eq. ( p.3| ). In the SM we obtain the 
following non-resonant branching ratios, denoted by Bnr, {£ = e,/i): 



Bnr{B 

BnriB - 

Bnr{B - 

Bnr{B - 

BnriB — 



-^ Ktt 

^ K*e+e 

K*T+T- 



5.7-10"^ 
1.3-10-^ 
2.3 • 10"^ 
1.9- 10"^ 
1.9-10"^ 



ABr,.r 



k27 



f^,±6,+^,±l,±2)%, 

+7 ~ 

=■=' 5-3 '-0.2 ? 



r+22 ^7+4 +o-\±i)% 



Ai3™. = (lf„+^,+12,+t,±3)% 



A g _ (+26 ,a +6 -0.7 ,91 
u^LJnr — 1-17) -'-U)-4 )+0.4 ' -"-^Z 



AB„ 



r+4 _A_/\ +13 +0.6 _|_o\Cr/ 



(5.1) 
(5.2) 
(5.3) 
(5.4) 
(5.5) 



The first error in the ABnr consists of hadronic uncertainties from the form factors. The other 
four errors given in the parentheses are due to the variations of m^, /x, mb^poie and as{mz), in 
order of appearance. In addition, there is an error of ±2.5% from the lifetimes tb |[T^. The 
scale-dependence of the branching ratio Bnr{B -^ K*e^e~) gives +12% and +1.4%, as /x is 
varied from /i = rrib^poie to /i = rrib^poiel'^ and /i = 2mb,poie, respectively, and we have taken 
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mw 


80.41 GeV 


mz 


91.1867 GeV 


sin^ 6w 


0.2233 


rric 


1.4 GeV 


'^bpole 


4.8 ± 0.2 GeV 


rrit 


173.8 ± 5.0 GeV 




+ 'ni,,poie 


^^ 


mb,pole_rrn,,^^lj2 


A (5) 


0.220i°:J]6l GeV 


a-1 


129 


asimz) 


0.119 ±0.0058 


VZVtb 


0.0385 


VtWt, / K.I 


1 



Table 6: Default values of the input parameters and the ±1 a errors on the sensitive parameters 
used in our numerical calculations. 



the larger of the two errors in this case to estimate the scale-dependence of this branching 
ratio. The largest parametric errors are from the uncertainties of the scale /z and the top 
quark mass, rrit. The large scale-dependence of the branching ratios reflects essentially that 
of the effective coefficients. To remedy this, one has to calculate the virtual corrections to the 
matrix elements of the partonic decays h —>■ si~^£~ to obtain perturbatively improved effective 
coefficients which are both scale- and scheme-independent [|7T]]. The exclusive decay form 
factors, obtained in the LCSR method including the radiative corrections, depend also on mj,, a^ 
and the renormalization scale /i. However, the various dependencies of the form factors are 
inadequate to compensate for the corresponding dependencies in the effective coefficients being 
used. We present in Fig. |^ the exclusive branching ratios calculated in the LCSR approach, 
obtained by adding the stated errors in quadrature. We also give, for the sake of completeness, 
the branching ratios for the inclusive decays B -^ {Xs,Xa)i^i^ . In calculating the theoretical 
dispersion on B ^ Xdi^i~ , we have varied the CKM parameters in the allowed range obtained 
from the CKM unitarity fits |jl2|. We have also listed the present experimental bounds on the 
exclusive decays B -^ {K,K*)e^e~ and B -^ {K , K*) fi~^ n~ , obtained by the CDF |0 and 



CLEO 1^ collaborations. Experimental upper limits on the inclusive decays B -^ X^e^e" 
and B -^ X^^^fi' are from the CLEO collaboration [33|. All experimental limits are 90% 
C.L., and for the sake of this figure we have averaged the branching ratios for the charged 
and neutral 5-meson decays, as the differences in their branching ratios are expected to be 
minimal theoretically. 

Figure D shows that the exclusive decays B -^ K*fi'^fi'' and B -^ K*e^e^ provide at present 
the most stringent bounds on the effective coefficients. While none of the experimental bounds 
has reached the SM-sensitivity, they do provide interesting upper limits on the parameter 
space of models with physics beyond the SM. We will discuss this point in detail below in 
the context of the SUSY models we are studying in this paper. We have also estimated 
the present theoretical precision on the quantity sq (zero of the FB-asymmetry) in the decay 
B -^ K*i'^i~' for i^ = e^ and i^ = fi^. Note, that due to the kinematics, there is no 
zero for the FB-asymmetry for the case B —* K*t^t^ . Theoretical uncertainties from the 
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Figure 3: Theoretical expectations for the exclusive decay branching ratios B{B — > K*i'^i~), 
B{B — >■ Ki~^i~), i^ = e^,fi^,T^, calculated in the LCSR method in the SM. For the sake 
of completeness, we also give the branching ratios for the inclusive decays B -^ Xgd^d' and 
B -^ Xdi^i~ in the SM, including the CKM dependence of the latter. Experimental upper 
limits (at 90% C.L.) are also shown: solid squares are from the charged B^ decays (and charge 
conjugate), circles from the decays of B^ (and charge conjugate) , and the empty squares are 
from the inclusive decays, averaged over the charged and neutral B decays. All experimental 
limits are from the CLEO ^d^ , [3^31/ and CDF / P^/ collaborations. 
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form factors and the four parameters discussed above, rrit, fi, mb,poie and as{mz), on the sq 
are estimated as: ±1%, ±0.3%, +14%/ — 7%, ±6%, ±4%, respectively. As discussed above, 
the form factor-dependent uncertainty for this quantity is minimal, and the main sources of 
errors are /i and rrib^poie- The reason of the marked scale dependence is again the lack of 
compensating perturbative corrections, in the absence of which the scale-dependence of the 
Wilson coefficients reflects itself in rendering sq rather imprecise. Adding the stated errors in 
quadrature, we estimate in the SM (fixing nib while varying fj, and ^^(m^)): 

So = 2.88+°:^^GeV' . (5.6) 

The actual dilepton mass distributions and the FB-asymmetry for the decays of interest in the 
SM will be given later, together with the corresponding estimates in some variants of SUSY. 

6 The decays B -^ {K, K*)£+£- in SUSY Models 

First studies of rare i?-decays B — * X^'-f and B -^ Xsi~^i~ in the context of MSSM were carried 
out in [jl6|, |17|, 0.Q Since then, these studies have been updated by taking into account progress 



in theory and experiments. We employ the following models to study the rare B -^ K^*H^ 
decays: (i) Minimal supergravity (mSUGRA), (ii) Relaxed SUGRA (rSUGRA), obtained from 

fni) 



mSUGRA by relaxing the universal scalar mass condition at the GUT scale [^, ^ ^ 
Minimal fiavor violating supersymmetric model (MFV) (in the sense that the fiavor violation 
is solely due to the standard CKM mechanism and resides in the charged current sector) pT| , 
and (iv) the Mass insertion approximation (MIA) ||2^. The last of these models serves as a 
generic supersymmetric extension of the SM having non-CKM fiavor violations. We do not 
consider models with broken R-parity and assume that there are no new phases from new 
physics beyond the SM, or, equivalently, that the constraints from the electric dipole moments 
of the neutron and charged lepton and indirect constraints from the decay B -^ X^j as well 
as other FCNC processes render these phases innocuous. This covers an important part of 
the supersymmetric parameter space, but not all. The issue of supersymmetric phases having 
measurable consequences in CP- violations in B and K decays and EDMs of the neutron and 
charged lepton is still far from being settled. As we have not studied CP-asymmetries in the 
decays B -^ {K, K*)i~^i~ , the neglect of additional CP-phases is not crucial to the analysis of 
the decay rates being presented here. 

The strongest constraint on the MSSM parameter space is coming from data on B ^ Xg'y 
[^, given in Eq. ( |1.1|) . In terms of the Wilson coefficients, this puts a bound on the modulus 
of Cy^ , given in Eq. ( |1.2| ) in the NLO approximation. The SM-based estimate of Cy*^ in the 
NLO precision is well within this range, which then restricts the otherwise allowed parameter 
space in the supersymmetric models we are considering. To be consistent with the precision 
of other contributions in B —>■ Xsi~^i~ , and for comparison with the rates and distributions 



in the SM, we work with Cj'^ i'mb,poie) in the LLA accuracy. This yields the bounds (at 95% 







C.L.) 

0.249 < |C7^^'^^^(/i = 4.8 GeV)| < 0.374 . (6.1) 



^There is a wrong sign in the chargino and ncutralino box matching condition in 1 16 . This sign discrepancy 
between jl^] and [18| has aheady been mentioned by the latter. We are grateful to T. Goto and F. Kriiger for 
clarifying this point. 



We remind at the outset that the theoretical uncertainties in the decay rates are estimated 
by us to be typically ±35%. Hence SUSY-searches in S — > {K, K*)£~^£~ will be unambiguous 
only for drastic SUSY effects. 

To illustrate generic SUSY effects in B ^ (K, K*)£^£~ , we start by assuming ICj'^^] ~ 
ICrf^l allowing for two possible solutions, ^7^^ < (SM-like) and Cj^^ > (allowed in SUSY 
models). We also fix the other two coefficients Cg and Cio to their respective SM values. We 
show the dilepton invariant mass distributions for B — >■ Kfi~^fi~ and B — >■ K*fi^fi^ decays in 
Figs. ^ and ^(a), respectively. The FB-asymmetry for B -^ K*fi~^fi~ is shown in Fig. ^(b). 
These figures present a comparative study of the SM- and SUSY-based distributions, and the 
attendant theoretical uncertainties associated with the long-distance effects. For the latter, we 
have used the KS-approach ||2^ and the AMM-approach [Q to estimate the resonance-related 



uncertainties. These figures illustrate that despite non-perturbative uncertainties, it will be 
possible to distinguish between the SM and a theoretical scenario in which the magnitude of 
the effective coefficients are similar, but C^'^^ has the "wrong sign". For the dilepton invariant 
mass, this reverses the sign of the interference term involving 'Re{{Cj'^^)* ■ Cg'^^) which leads to 
significant difference in the decays B -^ K*£^£~ . More striking deviation from the SM predic- 
tion is found in ApB for the models in which the condition Eq. ( 4.35 ) is not satisfied, resulting 



in a FB-asymmetry which remains negative below the J/?/'-resonance region. This would be a 
drastic deviation from the SM, which can not be fudged away due to non-perturbative effects. 
Interestingly, the situation Cy^ ~ — Cy^^^- is met in a number of SUSY models as discussed be- 
low. In addition, in a general flavor- violating supersymmetric model, also the other two Wilson 
coefficients (Cg and Ciq) may have either sign. In this case, the FB-asymmetry in i? — >■ K*£~^£~ 
may have a functional dependence on the dilepton mass which is characteristically different 
than the ones obtaining in the SM and SUGRA models, as shown below. 

More elaborate changes from new physics (NP) in the values of the relevant Wilson coeffi- 
cients can be taken into account by the (correlated) ratios, (i = 7, 9, 10): 






^iW — nSM ~ pSM ' y^-^) 



which depend on the renormalization scale (except for Cio), for which we shall always take 
/x = rrih^poie- The experimental constraint from B -^ Xg'y given in Eq. ( |1.1[ ) now translates into 
the bound 

0.80< |i?7(/i = 4.8 GeV)| < 1.20 , (6.3) 

where the coefficients are understood to be calculated in the LLA precision. In the numerical 
estimates, we have used Bgi = (10.4 ± 0.4)% for the average semileptonic branching ratio, and 
have set the heavy quark expansion parameters to the values Ai = —0.20 GeV^ and A2 = 0.12 
GeV^. The allowed values of the other two ratios Rg and Riq are taken from the literature for 
the mSUGRA and rSUGRA models |2|, ||], and for the other two models, MFV and MIA, we 
have calculated them. In particular, in the MIA approach, large enhancements are anticipated 



in the branching ratio B{B ^ Xs£ £~) in some allowed region of the parameter space 123 
These enhancements, suitably modified by the form factors, are also present in the branching 
ratios for the exclusive decays B — >■ {K, K*)£^£~ . However, as shown in Fig. |^, some of these 
branching ratios are bounded quite stringently, in particular, for the decays B —^ K*e^e~ and 



B -^ K* jj,^ jj, [^, ^. Assuming Rj in the allowed range, we shall work out the constraints 
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Figure 4: The dilepton invariant mass distribution in B —>■ Kfi~^fi^ decays, using the form 
factors from LCSR as a function of s. Solid curve: SD + LD using Ref. f^, dashed curve: 
SD + LD using Ref. /|^/; dotted: pure SD; long-short dashed curve: SD + LD using Ref. ^ 
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Figure 5: The dilepton invariant mass distribution (a), and the normalized FB-asymmetry (b) 
in B ^ K*fi^fi^ decays, using the form factors in LCSR as a function of s in the SM. Solid 
curves: SD + LD according to Ref. j&W, dashed curves: SD + LD using Ref. WW; dotted: pure 



SD; long-short dashed curves: SD + LD using Ref. MW with €■ 
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on the effective coefficients Cg and Cio (equivalently Rg and -Rio). Based on this analysis, we 
shall show the dilepton invariant mass spectra and the FB-asymmetry in some representative 

cases. 

6.1 B -^ {K, K*)£+i- in SUGRA models 

We shall consider here both the minimal and restricted SUGRA models (mSUGRA, rSUGRA). 
The parameter space of these models may be decomposed into two qualitatively different 
regions, which can be characterized by tan/3 values. For small tan/?, say tan/5 ~ 2, the 
sign of Cf'^^ is the same as in the SM. Here, no spectacular deviations from the SM can 
be expected in the decays B — > {K,K*)i^i~ . Given the theoretical uncertainties shown 
earlier by us, we think that it would be very difficult to disentangle any SUSY effects for this 
scenario in these decays. For large tan/5, the situation is more interesting due to correlations 
involving the branching ratio for B —>■ Xgj, the mass of the lightest CP-even Higgs boson, rrih, 
and sign{fisusy), appearing in the Higgs superpotential. In this case, there are two branches 
for the solutions for ruh and B{B —>■ Xsj). The interesting scenario for SUSY searches in 
B -^ {K, K*)i^i~ is the one in which sign{ij,susy) and nih admit Cj'^^ to be positive. For 
example, this happens for tan/3 > 10, in which case ruh = (115-125) GeV and Cj° is positive 



and obeys the B -^ Xg'j bounds |^. Following the generic case shown earlier, one expects 
a constructive interference of the terms depending on Cj^ and Cg in the dilepton invariant 
mass spectra. For the sake of illustration, we use 

i?7 = -1.2, Rg = 1.03, Rio = 1.0 , (6.4) 



obtained for tan/? = 30 pOf, as a representative large-tan/? solution, to study the effects on 
our observables. We find that in the low-g^ region the branching ratio for B — > Kfi^^~ is 
enhanced by about 30% compared to the SM one, as shown in Fig. p. This enhancement 
is difficult to disentangle from the non-perturbative uncertainties attendant with the SM- 
distributions (shown as the shaded band in this figure). The dilepton mass distribution for 
B — > K*n~^n~ is more promising, as in this case the enhancement is around 100%, see Fig. |^, 
and this is distinguishable from the SM-related theoretical uncertainties (shown as the shaded 
band in this figure). Note that the resulting branching ratios are consistent with the present 
experimental upper bounds on these decays given earlier. The supersymmetric effects presented 
here are very similar to the ones worked out for the inclusive decays B — >■ Xsi^£~ ||2^, where 



enhancements of (50-100)% were predicted in the low-g branching ratios. The effect of -Ry 
being negative is striking in the FB asymmetry as shown in Fig. |], in which the two SUGRA 
curves are plotted using Eq. ( |6.4D (for Rj < 0) and by flipping the sign of Rj but keeping 
the magnitudes of Ri to their values given in this equation. Summarizing for the SUGRA 
theories, large tan/? solutions lead to Cj^ being positive, which implies that FB-asymmetry 
below the J/'j/'-resonant region remains negative (hence, no zero in the FB-asymmetry in this 
region) and one expects an enhancement up to a factor two in the dilepton mass distribution 
in B ^ K*e+e- and B -^ K* ^+ ^-^ . 
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Figure 6: The dilepton invariant mass distribution in B —>■ K^i^ ix~ decays, using the form 
factors from LCSR as a function of s. All resonant cc states are parametrized as in Ref. j^/ . 
The solid line represents the SM and the shaded area depicts the form factor-related uncer- 
tainties. The dotted line corresponds to the SUGRA model with Rj = —1.2, Rg = 1.03 and 
Rio = 1. The long-short dashed lines correspond to an allowed point in the parameter space of 
the MIA-SUSY model, given by R-j = —0.83, Rg = 0.92 and Riq = 1.61. The corresponding 
pure SD spectra are shown in the lower part of the plot. 
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Figure 7: The dilepton invariant mass distribution in B ^ K*fi^fi~ decays, using the form 
factors from LCSR as a function of s. All resonant cc states are parametrized as in Ref. W^ 
The legends are the same as in Fig. [^. 
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Figure 8: The normalized forward-backward asymmetry in B ^ K*^^fi'' decay as a function 
of s, using the form factors from the LCSR approach. All resonant cc states are parametrized 
as in Ref. j^/ . The solid line denotes the SM prediction. The dotted (long-short dashed) lines 
correspond to the SUGRA (the MIA-SUSY) model, using the parameters given in Eq. \6.4( ) 
(Eq. (\6. 1 (\ )) with the upper and lower curves representing the C-{^ < and C^^ > case, 
respectively. The dashed curves indicating a positive asymmetry for large s correspond to the 
MIA-SUSY models using the parameters given in Eq. ( 6.1]\ ), i.e. the "best depression scenario" 
with Cio > 0. 



6.2 B -^ (K, K*)i+i- in MFV-SUSY Model 



The MFV-SUSY model is based on the assumption of minimal flavor violation. Here, quarks 
and squarks are aligned so there is no flavor-changing q — q' — {Z, 7, g) vertex and the charged 
one, d — u — x^, is governed by the CKM matrix. As a consequence, in this model neutralino- 
down-squark and gluino-down-squark graphs do not contribute to either 6 — >■ 57 or 6 — >■ si^i~ 
transitions. In addition to the charged Higgs-top graphs, chargino-up type squarks loops 
with a light stop ti, and the ly^-top quark loops, present in the SM, give the dominant 
contribution. While not holding generally, the assumptions in the MFV-SUSY model are 
valid over an important part of the minimal supersymmetric parameter space |^ . They have 
the simplifying feature that the dominant supersymmetric effects remain confined to charged 
current transitions and relatively easy to test experimentally due to well-defined correlations 
in several measurable quantities involving FCNC transitions [^ |12]. 

As is well-known []72[, in the two-Higgs doublet model of type II (2HDM model II), which is 
embedded also in the MFV-SUSY construct, the charged-Higgs contribution is always additive 
to the SM, i.e., C7^''^(2HDM) < 0, yielding a lower bound on the charged Higgs mass mH± 



(almost) independent of tan/?, above tan/3 > 1. In MFV, the x^ 
if ^ — t contribution, with a large positive contribution to Cj'^ . 



- ti loop can compensate the 
We scan over the parameter 
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space in the range 55 GeV < mH± < 1 TeV, < M2, \^susy\ < 500 GeV, where fisusy is the 
bihnear Higgs couphng in the superpotential and M2 is the gaugino soft breaking mass. We 
use rriq = m^^ = 1 TeV, where rriq denotes the (degenerate) masses of other than top squarks, 
and fix rric, = 50 GeV to its lower bound. We reject too hght charginos, demanding 'm^± > 70 
GeV, and also solutions which do not satisfy the bounds from the branching ratio on B ^ Xg'-f. 
The chargino contribution to Cy'^^ decreases for larger values of m^^ and we therefore keep it 
to its minimal value m^^ = 70 GeV. We have chosen a stop mixing angle O^ = ±27r/5, i.e. the 
light stop ti = cos OftL + sin 9^iji is almost right handed. 

For small tan/?, for which we again take tan/5 = 2, we find that the ratio R7 remains 
positive, i.e. Cj^ < 0, and lies within the experimentally allowed bounds from B -^ Xg'j, 
and the other two ratios are in the range 0.98 < Rg < 1.07 and 0.79 < Rio < 1.15. For 
large tan/3, taken to be 30, just as in the SUGRA models discussed earlier, Cj'^^ changes sign 
(-R7 < 0). The ratios -Rg and -Rio are again always positive but now Rg is almost identical to 
1, and RiQ tends to lie below the SM-value. Numerically, we find the ranges 0.99 < -Rg < 1, 
0.93 < -Rio < 1.02. The maximal (minimal) value of tan/3 found for -R7 > 0(< 0) is 5 (20). In 
contrast, a no mixing choice 9^ = ±7r/2 or ti ~ iji yields, for both tan/3 = 2 and 30, C^'^^ < 
or equivalent ly -R7 > 0. 

In general, in the MFV model, SUSY effects on Cg and Cio are much smaller than the 
corresponding one on Cj'^^. A large value of tan/3 helps Cj*^^ to satisfy the B — > Xs7 bounds 
but admits a sign opposite to the one in the SM. Dominant SUSY contributions to Cg and 
Cio are due to the charged Higgs exchange and are suppressed as ~ l/tan^/3, for large tan/?. 
Chargino effects in Cg and Cio increase for larger values of the ratio M2/\fJ,susy\ > 1- Using 
the central values of the parameters and the LCSR form factors, the maximal non-resonant 
branching ratios in the MFV are found for the ratios -R7 = —1.2, -Rg = 1.0 and -Rio = 1.02: 
jSrnaxMFV^Q _^ K ^+ ^-) = 7.5 X 10"^ and BZ^'^^^^i.B -^ K*n+i2-) = 3.2 x 10"^ While 
larger than the corresponding branching ratios in the SM, they are compatible with the present 
experimental bounds |34, 3^. Our findings in the MFV-SUSY model are very similar to the 



SUGRA case and in agreement with |T^ for the inclusive decays. As the values of Ri for 
the maximal non-resonant branching ratios in the MFV model are almost identical to their 
SUGRA- model counterparts given in Eq. ( |6.4|) , for which we have shown the dilepton invariant 
mass spectra and FB-asymmetry, we refrain from showing the corresponding figures for the 
MFV case. 

6.3 B -^ {K, K*)£+i- in the MIA Approach 

The minimal insertion approach aims at including all possible squark mixing effects in a model 
independent way. Choosing a. q,q basis where the q — q — x'^ and q — q — g couplings are flavor 
diagonal, flavor changes are incorporated by a non-diagonal mass insertion in the q propagator, 
which can be parametrized as {A,B =Left, Right) [l73| 



/ <rup,down\ \''''ij ) A,B /^ r\ 

{Sif )a,b = — 2 ' (6-5) 



m-g 



where (m"]'' "^"')a,b ^^^ the off-diagonal elements of the up (down) squark mass squared matri- 
ces that mix flavor i and j, for both the right- and left-handed scalars, and m? is the average 
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squark mass squared. The sfermion propagators are expanded in terms of the Ss. The Wilson 
coefficients have the following structure {k = 7, 9, 10): 

C, = d'' + Ct' + Cf ^^, (6.6) 



where C^^"^ is given in terms of (5"f' "^"')'ab ^P ^*^ ^*o mass insertions |^, and Cf^^""^ being 
the SUSY contribution in the basis where only flavor-diagonal contributions are allowed. It 
is tacitly assumed that the 6s are small and this defines the theoretical consistency of this 
approach which has to be checked a posteriori. 

The MIA-SUSY approach has been recently used in the analysis of the decays B — > Xsi~^i~ 
[^, taking into account the present bounds on the coefficient Cj'^ (rriB) following from the 
decay B -^ Xg'y. The other two coefficients C^^^ and C^^^^ are calculated by scanning over 



the allowed supersymmetric parameter space ||23[. For ^susy — —160 GeV, nig ~ m^ ~ 250 
GeV, mi^ = 90 GeV, m^ ~ 50 GeV, these coefficients are expressed as: 

C^^ = 1.75{61^,)ll - 8.25{6:^,)lr • (6.7) 

Of these, the mass insertions ((^23)^^ ^^^ (^23)^^ ^^^ related by a CKM rotation and the 
bound on one implies a similar bound on the other. One may have marked enhancement or 
depletion in the branching ratios for the decay B -^ Xsi~^i~. Note also the large numerical 
coefficient of (^23)^^ i^ ^^^ expression for Cf^^"^. For the parameters for which Eq. (|6.71 ) holds, 
the diagonal-SUSY contributions to Cg and Cio are: Cc>'^''"^{mB) = -0.35, Cw'^'"^ = -0.27. 
Depending on the value of (523)lh and (523)ll, the coefficient C^q^^ may easily overcome the 
SM- and the diagonal-MSSM-contributions in this coefficient, changing the overall sign of the 
FB-asymmetry. This feature is a marked difference between this scenario and the competing 
ones, namely SUGRA and MFV, where Cio remains close to the SM value (see Table 1). This 
feature has been noted already in [E31 in the context of the FB-asymmetry in the inclusive 



decay B -^ XJ+£-. 

To maximize the effects in this general flavor- violating supersymmetric context, several 



special cases have been studied in Ref. |^ in detail. We shall discuss the following three 
scenarios from this work:0 

1. "Best enhancement scenario" for the branching ratio B{B -^ Xsi^i~), which corresponds 
to the choice Cj"^ = 0.445, (6^3)11 = (^23)^^ = -0-5 and ((5^3)Li? = 0.9; 



2. "Best enhancement scenario with Cj^ < 0", corresponding to using Cj^ = —0.445 
iSl3)LL = -0.5, i61^3)LL = -0.1 and i6^3)LR = 0.9; 



3. "Best depression scenario", corresponding to Cj'^^ = —0.25, ((^23)^^ — 0-^' ("^23)^^ — 0-1 
and (^2^3)^^ = -0.6. 



^The specific values given above for the mass insertion parameter {S23)ll have been kindly provided to us 
by Ignazio Scimemi. We also draw attention to several misprints in the Tables given in |23] and trust that an 
Erratum is being issued by the authors of Ref. p3] . 
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With these choices, drastic effects in the branching ratios and the FB-asynimetry have been 
predicted for the decays B — >■ Xsi~^i~ , as displayed in Figs. 5-8 in Ref. [^. To wit, in the first 



scenario listed above, enhancements as large as a factor 5 are admissible in B{B -^ Xse~^e~) 
and even higher, 6.5, in B ^ Xsfi~^fi~ . 

We shall largely follow this analysis here in discussing the decay characteristics of the 
exclusive decays B — > {K, K*)i~^i~ but would like to add a dissenting remark concerning the 
coefficient Cj'^^{mB)- We recall that the extremal values used for Cj^^{mB) in p3[| correspond 
to using the 99% C.L. limits on B{B -^ X,7), which give the bounds 0.252 < jCy^^j < 0.445 in 
the NLO approximation. This procedure allows a much larger range for the ratio Rj than the 
one given in Eq. ( |6.3|) , which is then partly reflected in the branching ratios for B -^ Xsi~^i~. 

We argue that even with this more restricted range of Cy^ , the two "Best enhancement 



scenarios for B -^ Xsi^i~" of Ref. ||23| alluded to above give too large branching ratios for 
the exclusive decays being studied here. To be specific, in the first scenario, the parame- 
ters given above translate into Rg = 1.26 and -Rio = 2.84.|^ The central values of the form 
factors calculated here in the LCSR approach then lead to the following branching ratio: 
j^rnax,MiA^j^ — >■ K*fi~^fj,~) = 11.5 X 10~^, which is approximately 3 times larger than the recent 



CDF (90% C.L.) upper limit on this quantity ||], 



S(5° -^ i^*W") < 4.0 X 10"^ . (6.8) 

The B ^ K transition in this scenario is likewise enhanced, yielding a branching ratio 
j^'^^^MiAi^j^ _^ i^/i+yU^) = 3.2 X 10^^, which is typically a factor 5 larger than the SM branching 
ratio, but still compatible with the experimental upper limit, B{B^ -^ K^fi~^fi~) < 5.2 ■ 10~^ 



JM. Hence, the present experimental upper bound on B ^ K* jj,'^ jj, provides non-trivial 



bounds on Cg and Cio, equivalently on R^ and i?io, which we now proceed to work out. 

6.4 Bounds on Cg and Cio from present data 

The branching ratios B -^ {K, K*)£^£~ can be expressed as quadratic equations in the co- 
efficients Cy*^^, Cg and Ciq. Given the branching ratios (equivalently upper bounds), these 
equations can be solved numerically and yield the allowed contours in the Cg-Cio plane. For 
working out the constraints, we use the experimental bound in Eq. (|6.8| ) and the following 
expression which follows from Eq. ( |4.27| ): 

B{B ^ K*f,+f,-) = a^;^:^\Cr''''\'+bP\Cg\'+ci:^\C,o\'+dPCr'''Cg+e'-r^^^ • 

(6.9) 
The coefficients a^*, •••, Qk* ^^^ tabulated in Table ^ using the central values of the B -^ K* 
form factors in Table ^ and the maximum and minimum values of the same given in Tables ^ and 
1^, respectively. Of these, the coefficients bj^V and Cj^V coincide if one neglects the £='=-masses. 
The superscript on these coefficients is a reminder that only non-resonant contributions are 
included. 

The quadratic equation in ( |6.9D is solved numerically for the two distinct situations Cj'^^ < 
(SM-like) and Cj^ > (new physics scenario) in the experimentally allowed range for Cy*' 



given in Eq. (|6.1| ). The resulting 90% C.L. allowed contours are shown in Fig. ^ and Fig. 10, 



We neglect the effect from the RG running from /i = ms (used in |E3[) to /i = nib^poie used by us. 
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Anr) 


Anr) 


„{nr) 


Anr) 


Anr) 


r.(nr) 
JK* 


(nr) 


FF (central) 

FF(max) 

FF(min) 


21.295 
28.183 
16.795 


0.502 
0.630 

0.417 


0.500 
0.633 
0.416 


3.530 

4.577 
2.864 


1.434 
1.859 
1.164 


0.413 
0.520 
0.343 


0.148 
0.183 
0.125 



Table 7: Coefficients of the non resonant branching ratio B{B -^ K*fi^fi~) in units of 10^''' 
in the decomposition as in Eq. ( \6.^) , integrated over the full q^ range for different sets of form 
factors given in Tabs. [^-[^. 

respectively. The solid curves in these figures are obtained by using the central values of the 
form factors and the inner and outer dashed curves represent the maximal and minimal allowed 
values of the same, respectively. Note that the loosest bounds emerge from the minimal allowed 
values of the form factors. Also, in working out the constraints shown in these figures, we have 
fixed \Cy'^ I = ICy^j^l = 0.249 in the allowed range given in Eq. (|6.1|) , as this gives for both 
the cases (Cy^ < and Cj^ > 0) the loosest bounds on Cg and Ciq. This can be seen 



in Figs. |Tl| and Fig. |T^ drawn for Cj'^^ 



eS 



< and Cj'^ > 0, respectively, where we show the 



dependence of the bounds in the Cg-Cio plane on the experimentally allowed range for IC7 



effi 



given in Eq. (|6.1|) . In these figures, we use the minimum values of the form factors given in 
Table ^for reasons given above. In Figs. ^and|TT|, we also show the SM-point (see Table 1) and 
the SUSY-MIA points for the "Best enhancement scenario with Cy'^^ < 0", corresponding to 
Cc,(mFi) = 5.0, Cin = —12.5, and the "Best depression scenario with Cf'^^ < 0", corresponding 

0.2, |]23|. We note that the "Best enhancement scenario with Cy'^^ < 0" 



to Cg^niB 



'5 ^10 ■ 

3.2, C 



10 



is ruled out by data. The other MIA-SUSY point, as well as the SM, are both well within the 
experimental bound. The SUSY-MIA point corresponding to the "Best enhancement scenario 
with Cy^ > 0" of Ref. [^ is shown in the C9-C10 plane in Figs. and |12|. This corresponds to 
the point Cg{mB) = 5.5, Cio = —13.2. As anticipated, this "Best enhancement scenario with 
Cj'^ > 0" is convincingly ruled out by the experimental upper bound on B{B -^ K*ix^fj.~). 
The analysis shown in Figs. ^0 holds for all models discussed here in this paper in which 
the SD-physics can be encoded in terms of the three real Wilson coefficients C^^^ , Cg and Cio- 
The point we wish to stress is that existing data on S — >■ K*ij,~^ii~, in conjunction with the 
branching ratio B{B -^ Xg'j), provides non-trivial constraints on Cg and Cio- 

Illustrative examples of the dilepton invariant mass spectrum in the decays B -^ K^'^^~ 
and B -^ K*^^fi~ in the MIA approach are shown in Figs. |^ and |^, respectively. They have 
been calculated for the following values: 



R7 = ±0.83, Rg = 0.92, R 



10 



1.61 



(6.10) 



which are allowed by the present experimental bounds. The characteristic difference in this 
case, as compared to the SUGRA and MFV-SUSY models, lies in the significantly enhanced 
value of CiQ. 



As already mentioned, a characteristic of the MIA approach is that the sign of Ciq {C 



SM 
10 



< 



0) depends on the quantities (^23)^^ ^^d {S23,)ll- In particular, the large number in front of 
(^23)l-R in C'lo, obtained for the specific values of the SUSY parameter space, could change the 
sign of this Wilson coefficient. This has no effect on the dilepton invariant mass distributions, 
as they depend quadratically on Cio, but it would change the sign of ApB in 5 — > K*i~^£~. 
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Figure 9: Bounds on the coefficients Cg{mB) and Ciq resulting from the experimental upper 
bound B{B° -^ K*^i2+fi-) < 4.0 x 10"^ (at 90% C.L.) §^ and Cj^'^ifi = 4.8 GeV) = -0.249 
from the bounds given in Eq. ( |g. j| / The SM-point and two representative points in the SUSY- 
MIA approach from Ref. [2^J are also shown. The three curves correspond to using the central 
values of the form factors (solid curve), the minimum (outer dashed curve) and maximum 
(inner dashed curve) allowed values discussed in Sec. 3. 
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Figure 10: The same as Fig. ^ but for the solution with C-j^ = 0.249. The point MIA^est 
corresponds to the "best enhancement scenario" of Ref. WW, discussed in the text. 
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Figure 11: The same as Fig. ^ hut showing the dependence of the hounds on the experimentally 
allowed range for \Cj^ \, 0.249 < \C'j^ \ < 0.374, with the form factors fixed to their minimum 
values given in Tahle^. 
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Figure 12: The same as Fig. UB, hut with Cj > 0. 
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a 


h 


c 


d 


e 


/ 


9 


h 


J 


k 


K 
K* 


0.193 
13.119 


0.068 
0.197 


0.068 
0.196 


0.230 
1.760 


0.163 
0.995 


0.097 
0.236 


0.045 
0.083 


0.943 


0.089 


0.061 



Table 8: Coefficients in units of 10 '^ defined in Eqs. ( \6.13i ) and ( ^. 14 ) in the KS prescription 



To illustrate this, we use the parameters close to the so-called "Best depression" scenario p3| , 
corresponding to the following values 



R^ = ±0.83, Rg = 0.79, R 



10 



-0.38 



(6.11) 



and plot the resulting normalized FB asymmetry in Fig. ^. The positive FB-asymmetry in 
B -^ K*i^i~ (as well as in i? — >■ Xsi^i~ shown in ||23|) for the dilepton invariant mass below 
the resonant J/ip region is rather unique, as none of the other models considered here (SM, 
SUGRA and MFV) admit solutions with positive Ciq. 

Finally, to facilitate a model independent determination of the coefficients Cj*^ , Cg, and 
Cio from the decays B -^ {K,K*)i~^i~, we write down a parametrization of the partially 
integrated branching ratios and FBA in the low s region. Using, for the sake of definiteness. 



0.25 GeV^ Sr 



8.0 GeV , the partial branching ratios ABx and the corresponding 



FB-asymmetry AAfb can be expressed as (X = K, K* 



AB 



X 



ds 



dB{B -^ Xfi+fj- 
ds 



axlCj^ I +&x|C*9| +cx|C'io| 



dxC-j^ Cg + exCj^ 



(6.12) 



fxCg + gx (6.13) 



AAfb = tb 



ds 



dApBJB -^ K' 
d7~ 



^+^ ) 



Cio{hxC-j^ 



jxCg + k 



X, 



(6.14) 



Numerical values of the coefficients are given in Table § They have been obtained by using 
the central values of the form factors and other parameters given in Table 13 and Table |^, 
respectively. Specifying a model by the effective coefficients C7'^^{mB), Cg'^ (mB) and Cio 
enables one to obtain readily the predictions for ABx and A^ in this model. In the SM, 
we estimate ABk = 2.90 ■ 10"^ ABk- = 7.67 ■ 10"^ and A^^.^ = -0.71 ■ 10"^ yielding 
AAfb = AAfb/ ABk* = -9.2%. The branching ratios for the decays B -> {K, K*)e+e- 
are practically identical. Typical theoretical errors on these quantities, obtained by varying 
the form factors and the parameters rrit, nib, f^ and Aqcd in the ranges discussed earlier and 
adding the individual errors in quadrature are ±30% for AB and ±38% for A^. However, the 
branching ratios and the FB-asymmetry may be significantly enhanced (or depressed) in some 
variants of the supersymmetric models discussed. With 0(10*^) BB events anticipated at the 
B- factories and HERA-B, and much higher yields at the Tevatron and LHC experiments, these 
rates and asymmetries will allow precision tests of the SM and may indicate the presence of 
new physics. 
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7 Summary and Concluding Remarks 

Before summarizing our results, we would like to comment on the contributions from the 
helicity-flipped counter-parts of the SM operators C7, O9 and Ciq: 

^7 = ^sa^^m.LbF^-' , (7.15) 

^9 = T^.snl'bnh/ , (7.16) 

^10 = T^snl'-bRh.l^i . (7.17) 

In an enlarged operator basis including these and the SM-operators, the various distributions 
for the decays of interest can be obtained from the substitutions Ci ^ Ci + C[ {i = 7, 9, 10) in 
the matrix elements and the auxiliary functions Eqs. ( [4.6|) -( ^rT7D for B —>■ K, and for B -^ K* 



in the terms which are proportional to the form factors V and Ti. In the remainder of the 
B -^ K* amplitude, the contribution of the helicity-flipped operators enters with the opposite 
sign, i.e., Ci ^ Q - C[. 

We note that in all models with minimal flavor violation, like the SM, 2HDM, and MFV, 
the contributions of the flipped operators Oj^ ^q vanish in the tUs -^ limit. In the general 
non-diagonal MSSM scenarios, there are finite contributions even for a vanishing s-quark mass 
due to the neutralino-gluino-down-squark loops. However, under the assumption that no large 
cancellations happen, we can conclude from the data on B{B -^ Xs'j) which bounds ICj*^^ + 
Cf^'l'^ that Cf^' must be small compared to Cj'^^. Further, neglecting box diagrams, the helicity 
structure of the (penguin)-loops responsible for Cg ^q can be related to the ones of the flipped 
photon penguin C^^' and hence is suppressed as well. We also note that we have neglected 
the effects of the neutral Higgs exchanges, which may lead to some inaccuracies for the decay 
B — > {K, K*)t^t~ in some parts of the SUSY parameter space. They are insignificant for the 
decays involving the {K, K*)fi~^fi~ and {K, K*)e^e~ states, where most of the experimental 
searches will be concentrated. 

We summarize our results: We have undertaken an improved calculation of the form fac- 
tors in the decays B — ^ {K, K*)i~^i~ in the light cone QCD sum rule approach. Using this 
framework, we have calculated the partial branching ratios, dilepton invariant mass spectra 
and the forward-backward asymmetry for these decays in the context of the SM. We have also 
undertaken a comparative study of the phenomenological profiles of these decays in a number 
of supersymmetric models. These include the SUGRA models, minimal-flavor- violation SUSY 
model, and a general flavor-violating SUSY framework using the mass insertion approxima- 
tion. The role of the forward-backward asymmetry in the decays B —>■ K*i^i~ in searching 
for new physics is emphasized. We show that the large- (tan /5) solution in the SUGRA models, 
but also some parameter space of the MIA model, yield FB-asymmetries, which are strikingly 
different from the SM. In particular, the value of the dilepton invariant mass for which the 
FB-asymmetry may become zero, sq, may provide a precision test of the SM. A simple analytic 
expression for sq is derived, and we have argued that the form factor dependence in sq cancels in 
the large energy expansion approximation. We have analyzed the present data on B ^ Xg'j and 
existing limits on the decays B — > {K, K*)i^i~ to put bounds on the coefficients Cg and Ciq. 
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While these bounds do not yet probe the SM, they do provide non-trivial constraints on exten- 
sions of the SM. In particular, the "Best enhancement SUSY-MIA scenarios "for the branching 
ratios B{B -^ Xsi^i~), shown for some chosen supersymmetric parameters in Ref. p3 



are 



ruled out by the existing upper limit on the exclusive branching ratio B{B^ -^ K*^^^ ^~) pl |. 
Finally, we show the dilepton mass spectra and the FB-asymmetry for illustrative values of the 
supersymmetric parameters and argue that the decays B -^ {K,K*)£~^i~ hold great promise 
in unraveling new physics. 
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